Vortex dynamics in the presence of a line of submicron holes along a superconducting 

microbridge 
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We measured and compared the electric field vs current density characteristics in the vortex state 
of two amorphous Nk>o.7Geo.3 microbridges, with and without a line of submicron holes patterned 
along the sample axis. The power dissipation in the perforated sample exhibits a crossover, being 
reduced at temperatures well below the superconducting transition temperature T c and unexpectedly 
enhanced close to T c . At low temperatures the holes are efficient artificial pinning centres and 
reduce the average vortex velocity. We argue that the dissipation enhancement close to T c is a 
consequence of a combination of the weakened pinning by the holes and an inhomogeneous driving- 
current distribution in their vicinity, which results in an increased average vortex velocity as well as 
in a channeling of the vortex motion through the holes. 

PACS numbers: 74.78.Na, 74.40.+k, 74.25.Qt 



The rapid development of nanostructuring methods in 
the last decade has enabled highly precise fabrication of 
artificial pinning centres (APCs) for vortices in super- 
conductors. It has already been proved that point- like 
inclusions in a superconducting film, such as perforating 
holes^ structural defects 3 - and magnetic dots^ can pin 
vortices. A regular two-dimensional array of APCs sta- 
bilises the vortex lattice against external driving forces, 
and results also in commensurability effects which scale 
with a matching magnetic field Bm — (f>o/S, where S is 
the area of the array unit cell and <f)Q the magnetic flux 
quantum. 1 At a magnetic field B > Bm vortices fill in 
the interstitial positions as well, where they are subject 
to intrinsic pinning originating from the structural de- 
fects. In a certain B range above Bm they may even 
form a commensurate multiple-flux-quanta lattice, as re- 
vealed from magnetisation measurements^ A less sym- 
metrical APC landscape, such as a rectangular array£ or 
a set of parallel continuous lines of a magnetic material, 6 
introduces anisotropies to the vortex transport. So far 
the effects of APCs have been found to be important at 
temperatures T rather close to T c , typically at reduced 
temperatures t = T/T c > 0.9, whereas at lower tempera- 
tures the intrinsic pinning dominated. Moreover, APCs 
invariably decreased the dissipation, apart from in a very 
recent report £ 

The simplest APC is a small perforating hole in a su- 
perconductor, having magnetic permeability larger than 
its diamagnetic surroundings. The most efficient pinning 
is obtained if the hole diameter is comparable to the 
magnetic-field penetration depth A rather than the co- 
herence length £, as one may expect at first sight. This 
was predicted by virtue of the Ginzburg-Landau (GL) 
theory for a single hole& and observed experimentally 
for arrays of holes^ The relevance of A for the artificial 
pinning potential is a consequence of the comparatively 
small supercurrent kinetic energy required for preserving 
the flux-quantisation condition around a hole not smaller 



than A. 

Local interactions of vortices with holes in the pres- 
ence of a driving current have not yet been studied in 
detail experimentally. These go beyond the cited match- 
ing phenomena and are important for understanding the 
function of holes in an arbitrary arrangement, which may 
have implications for possible design of future devices 
based on the manipulation of artificially pinned vortices. 
In this report we present such an investigation facili- 
tated by arranging the holes in a line centred along the 
main axis of a superconducting microbridge [lower inset 
to Fig. 1(a)], which eliminates the matching phenomena 
but still provides a sufficient signal arising from the pres- 
ence of the holes. In addition, we have minimised the 
influence of the intrinsic pinning by using the amorphous 
superconductor Nbo.7Geo.3j which is known for its very 
low background pinning>iS*ii 

The vortex motion was detected by recording the elec- 
tric field vs current density characteristics E(J) in three 
different temperature regimes, over a wide range of ap- 
plied currents, and B up to the upper critical magnetic 
field B C 2- Close to T c we observed a clear increase of 
the power dissipation by the holes over the whole range 
of B and J. This result can be explained plausibly by 
taking into account a local enhancement of the vortex 
driving force due to a spatial modulation of the driving- 
current density around the holes. To our knowledge such 
current-modulation effects have so far been disregarded 
in the interpretation of previous experiments carried out 
on perforated superconductors. Well below T c the holes 
pin the vortices and reduce the dissipation, which shows 
that their usefulness as APCs may extend to low temper- 
atures if the intrinsic pinning is weak. Our results imply 
that the interplay of hole pinning and inhomogeneous 
current drive determines whether the holes enhance or 
reduce the dissipation in the vortex transport. 

Two 210 fxva long, 5 fim wide and 20 nm thick mi- 
crobridges between two large contact pads were sput- 
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FIG. 1: E(J) of sample N (o) and sample P (•) for (a) T = 1.1 
K , B = 0.9 T, (b) T = 2.0 K, B = 0.6 T and (c) T = 2.5 
K, B = 0.25 T. The solid and dashed lines are plots of LO p/ 
expected at low and high t, respectively. Lower inset to (a): 
A sketch of sample P (not to scale) and the designation of the 
directions (see the text). Upper inset to (a) and insets to (b), 
(c): The initial parts of the E(J) in an expanded scale for a 
better view. 



tered onto the same Si/SiC>2 substrate as described in 
ReflToL Both samples had the same T c = 2.75 K and the 
same overall normal-state and superconducting proper- 
ties. Throughout the paper we use the nonperforated 
sample (sample N) as a reference for identifying the ef- 
fects of holes. Its transport properties were analysed in 
detail in RefUll and its GL parameters £(0) = 6.8 nm 
and A(0) = 1.15 pxs\ are taken as representative of both 
samples. The perforated sample (sample P) contained a 
line of equidistant holes along the central axis of the mi- 
crobridge (designated as x direction) as sketched in the 



lower inset to Fig. 1(a). The holes had a diameter of 
800 nm [slightly smaller than A(0)] and their centre-to- 
centre distance was 1.2 fim. A magnetic field was applied 
perpendicularly to the film plane and a dc current was 
passed in the x direction, so that the vortices traversed 
the sample in the y direction. The current sweep rate 
was 10 nAs -1 . The applied current is for sample P con- 
verted into the average current density Jp by calculating 
the (uniform) current density Jn for sample N and tak- 
ing the ratio of the normal-state resistances to determine 
Jp = 1.23Jat. We use this notation, together with En, 
Ep for the electric field in samples N and P, respectively, 
when referring to the E{ J) of the two samples specifi- 
cally. At high J the E(J) exhibit nonlinearities due to 
the flux-flow instabilities that are related to the nonequi- 
librium changes of vortex cores^i Here we concentrate on 
the close-to-equilibrium regime where these effects have 
not yet been developed and the vortex cores maintain 
their equilibrium properties. 

In Fig. 1 we plot typical E(J) representative of three 
temperature regimes: (a) T well below T c (1.1 K, t = 
0.4), (b) intermediate t (2.0 K, t = 0.7), (c) T close to T c 
(2.5 K, t = 0.91). The open circles show En(Jn), solid 
circles Ep(Jp), and the lines different theoretical pre- 
dictions of the Larkin-Ovchinnikov (LO) flux flow (FF) 
theory^ as discussed below. The scaled magnetic field 
b = B/B C 2 corresponds to 0.30 - 0.38, the B C 2 values 
being (a) 3.0 T, (b) 2.0 T and (c) 0.65 T, with the uncer- 
tainty of around 5 %. At low T [Fig. 1(a)] the E(J) of 
both samples reveal two different dynamic regimes: ther- 
mally activated magnetoresistance at J — * 0, followed at 
larger J by an E oc ( J — J c ) behaviour that implies FF 
against a background pinning potential with a depinning 
threshold J c ~ The solid line is a plot of the latter de- 
pendence, with the slope dE/dJ = p n b/0.9 equal to the 
low-i LO FF resistivity pf (p n is the normal-state resis- 
tivity), and J c chosen to obtain a fit to En(Jn)- The LO 
theory describes this part of E(J) reasonably well until 
the out-of-equilibrium nonlinearities in E(J) start to take 
place at large Jii As can be seen, the power dissipation 
in sample P is lower for ~ 10 % throughout the whole 
current-density range, suggesting an efficient pinning by 
the holes even far below T c . In Fig. 1(b) we show the 
E(J) at T = 2.0 K (t = 0.7), just at the boundary of the 
low-i and high-£ regimes, displaying a weaker reduction 
of Ep below En. Close to T c [Fig. 1(c)] the theoretical 
high-i LO pf (see Refllllfor details), shown by the dashed 
line, describes En(Jn) excellently starting from J = 
and up to the appearance of the nonlinearities mentioned 
before. However, over the entire current-density range 
the dissipation in sample P is larger than in sample N, 
unexpectedly and in contrast to the result at lower tem- 
peratures. The reason for this peculiar behaviour cannot 
be found in the pinning properties of holes. 

The magnetic field range over which En Ep is wide 
at all three characteristic temperatures. This is demon- 
strated in Fig. 2, where we plot the ratio p = Pp/Pn 
of the power density dissipated in samples P and N vs 
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FIG. 2: The ratio p = Pp/Pn of the power dissipated in 
samples P and N, integrated up to the appearance of the 
flux flow instabilities, vs b — B/B C 2- Sufficiently below T c 
the holes reduce the dissipation (p < 1) whereas close to T c 
the situation is opposite. 

b = B/B C 2- The integration Pn,p = J En,p<IJn,p is for 
each curve performed over the maximum region where 
the dissipation is not affected by the FF instabilities. 
We have checked whether p depends on the upper limit 
of integration, and we found minor numerical differences 
of the order of 10 - 15 % of the values shown in Fig. 2, 
with no change in the general shape of p(b). As can be 
anticipated from the E(J) shown in Figs. 1(a), 1(b), for 
t = 0.4 and t = 0.7 we find p < 1, i.e. the holes are active 
pinning sites and decrease the dissipation. Their relative 
contribution to the pinning becomes smaller as t and b 
increase, which is expected qualitatively and discussed 
later. Close to T c the dissipation in sample P is always 
larger, thus p > 1 with a maximum around b ~ 0.4. This 
enhancement implies a suppression of the artificial pin- 
ning by another effect which is fully manifested close to 
T c . 

In order to understand the above results one has to 
address the equation of motion for the vortices 

mu = F d - rju + F h , (1) 

where u is the vortex velocity, ry the vortex-motion vis- 
cosity, m the effective vortex mass per unit length^ 
Fd = 0oJ x z the driving force (z is the unit vector in 
the direction of the applied magnetic field), and F^ be- 
ing the repulsive force between the mobile vortices and 
those pinned by the holes. The spatial dependence of the 
z-th component of the vortex acceleration u is given by 
iti = u • Vu,. For our qualitative reasoning we have ne- 
glected other contributions (e.g. the interaction of mo- 
bile vortices with each other or with intrinsic pinning 
potential). Knowing the supercurrent distribution J(r) 
one can solve Eq. Q to find the vortex-velocity pro- 
file u(r). The generated electric field, the x component 
of which contributes to measured Ep, is calculated as 
E(r) = u(r) x B(r), where B(r) is related to the local 
vortex density n(r) = B(r)/<f>Q. When confining a pinned 
vortex each hole acts as a source of repulsion to other 
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FIG. 3: A qualitative consideration of the vortex trajectories 
(arrowed solid lines) around the holes at (a) T close to T c 
and (b) T well below T c . The modulated current is indicated 
by the arrowed dashed lines, the density of which represents 
the magnitude of J. Arrowed dotted lines in (b) depict the 
repulsion excerted by the vortices pinned by the holes to other 
vortices. The labels L and H denote regions of low and high 
vortex flow density, respectively. 



vortices, with pointing radially from the centre of the 
hole. The repulsive force between two vortices separated 
by r <C A is given by A _2 ln(A/r) and therefore weak- 
ens as A — > oo. Furthermore, the probability of a vortex 
being pinned by a hole becomes progressively smaller as 
A grows much larger than the hole diameter^ Hence it 
is reasonable to assume that decreases monotonically 
with increasing A and eventually becomes irrelevant in 
the limit of diverging A at T — > T c . In the following we 
discuss the results for t = 0.91 assuming Fft = 0. 

We note that the driving force F^ is not uniform 
around the holes because the supercurrent density J is 
spatially modulated. A schematic of the modulation of 
J, and consequently of Fj, around the holes is shown in 
Figs. 3(a) and 3(b) by the dashed lines. In the regions 
denoted by H the driving force is larger than far from the 
holes, while in the regions denoted by L the situation is 
the reverse. With F^ rj in the vicinity of T c , the vor- 
tices are accelerated parallel to F<; (perpendicular to J), 
thus the vortex trajectories accumulate in regions H as 
sketched in Fig. 3(a) by the solid lines. The second effect 
of the current modulation is that in regions H the vortices 
move faster than in the sample bulk and produce a larger 
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local electric field, while in regions L the opposite hap- 
pens. Although the holes themselves shorten the distance 
the fast vortices travel in regions H, a sufficient imbalance 
in favour of the number of these vortices, together with 
their large acceleration, may result in the total dissipa- 
tion in sample P being larger than in sample N. Therefore 
if a highly dense two-dimensional array of holes 1 * is in- 
tended to be used for enhancing the pinning of a sample 
there is no guarantee that this will work close to T c in the 
presence of a transport current, although the magnetisa- 
tion loops may widen up. Moreover, the current-induced 
vortex-velocity enhancement strengthens with increasing 
J, which sheds more light on the similar result of Ref0 
which also occurred at a relatively large applied current. 

The nonmonotonic b dependence of p > 1 may be 
linked to the competition between the enhanced current 
drive around the holes and the mutual repulsion of the 
mobile vortices, as explained below. The excess electric 
field relative to that far from the holes is estimated as 
AE ~ (f>oAuAn, where Ait and An are the effective ex- 
cess vortex velocity and density in regions H, respectively. 
While Au depends on J, n, m, and is independent of B, 
An also depends on the repulsive interaction between the 
mobile vortices, which is weak at low b and strong at high 
b. At low fraction b = B/B C 2 of the volume filled with 
the vortices they can all easily be channeled through re- 
gions H. Thus a larger vortex density results in the larger 
An and AE. On the other hand, An is at high b limited 
by the reduced space available for the channeling, which 
decreases AE and in turn results in a nonmonotonic p(b) . 

As temperature is lowered A becomes smaller and the 
holes start to pin more efficiently,- which explains quali- 
tatively the stronger reduction of the dissipation at lower 
temperatures (Fig. 2). The vortices trapped by the holes 
now repel the incoming vortices and change the vor- 
tex trajectories depicted in Fig. 3(a) to those shown in 



Fig. 3(b). The repulsive force Fh is indicated by the dot- 
ted arrows. In a simple model of pinning by the holes in 
dynamic conditions, a vortex remains pinned by a hole 
for some time until it is replaced by an incoming vortex. 
At low vortex density almost all vortices in the vicinity 
of holes are either pinned by them or scattered to pass 
through the regions L of low current density. Hence, the 
holes cause a noticeable reduction of the dissipation. As 
B increases the incoming vortices exert more force upon 
those pinned at the holes, reduce the pinning time, and 
the overall suppression of the dissipation decreases. This 
qualitative picture may explain the behaviour of p(b) at 
low temperatures. 

In conclusion, we have patterned a line of holes with 
a diameter close to A(0) along an amorphous Nbo.7Geo.3 
microbridge. A comparison of the measured E{ J) curves 
of the samples with and without perforation reveals an 
unusual crossover in the power dissipation close to and 
well below T c . Close to T c the artificial pinning is weak 
and an unexpected rise of the dissipation is observed. 
This is attributed to the inhomogeneous current distri- 
bution around the holes leading to a significant increase 
of the local vortex velocity. As temperature is lowered 
the pinning by the holes becomes stronger and eventu- 
ally suppresses the vortex velocity enhancement. In addi- 
tion, our weak background pinning has permitted, to our 
knowledge, the first observation the pinning properties of 
holes as artificial pinning centres far below T c . 
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